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Abstract 

The convergent iterative procedure for solving the groundstate Schroedinger equation is 
extended to derive the excitation energy and the wave function of the low-lying excited states. 
The method is applied to the one-dimensional quartic potential problem. The results show 
that the iterative solution converges rapidly when the coupling g is not too small. 
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1. Introduction 



Consider the iV-dimensional Schroedinger equation 

(-\V 2 + V(q)-E)ij(q) = (1.1) 

with 

q = (qi, 92, • • • , qw) (1-2) 

and 

v 2 = V — 

i=\ dqf 

Through a linear transformation of the coordinates qi, q^ qN, (11) 

can be applied to most of the nonrelativistic many body problems. Similarly, 
in the limit N — > oo, the same equation can also be extended to relativistic 
bosonic field theories. However, except in a few special cases, it is difficult 
to derive analytical solutions. Recently, for a class of Schroedinger equations 
we have succeeded in deriving a convergent iterative series solution for the 
groundstate[l-4]. These include the Sombrero shaped potential in any space 
dimension and arbitrary angular momentum[5]. In this paper, we discuss pos- 
sible extension to low lying excited states. 
Let 

MQ) = ^ S(q) (1-3) 
and E g d be the ground state wave function and the ground state energy. De- 
note the corresponding ones for the first excited state as 

^ ex = e~ s x and E = E gd + e. (1.4) 

We assume the potential V(q) to have high barriers and low valleys, so that the 
excitation energy e is very small. As an example, take q to be a one-dimensional 
space x and 

V(x) = 9 ^(x 2 - l) 2 (1.5) 
with g large. The first excited state would have an excitation energy 

e - e~^ 9 « 1. (1.6) 
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For clarity, we shall also assume that, in the general case of V(q), there is 
only one such low-lying state, as in the example of the one-dimensional quartic 
potential. 

Let the Schroedinger equations for ip g d = e~ s and ip ex = e~ s x be 

{-\v 2 + V-E gd )e- s = V (1.7) 

and 

(-^V 2 + V- E gd )e- s X = ee- s X . (1.8) 
Multiplying (1.7) by e~ s x and (1.8) by e~ s , we obtain from their difference, 

-iv • (e- 25 V X ) = ee~ 2S X . (1.9) 
Replace (1.9) by a series of iterative equations: for n > 1 

-iv • (e- 2S V Xn ) = e n e- 2S Xn -i (1.10) 

and when n = 0, 

e~ 5 xo = a trial function ^^a/ (1-H) 
which satisfies the orthogonality condition between ip g d and iptrial', i e. 

Je- 2S X od N q = 0. (1.12) 

When Xn-i is known, the iteration of (1.10) gives only the ratio Xn/^n- To 
perform the next order iteration it is necessary to separate e n and % n . For this 
purpose, we may choose a fixed point q° and set 

Xn(q°)=Xo(q°). (1.13) 

Clearly, different choices of g° will yield different sequences of x n (<?) and e n , 
as we shall discuss. 

For given e and x n _i, (1-10) can be viewed as an analog problem in 
electrostatics with a dielectric constant k and an external electric charge a n 
given by 

K = e- 2S and a n = e n e- 2S Xn -i (1-14) 
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Let E n be the corresponding electrostatic field and D n the displacement field. 
We write 

E n = -\v X n (1.15) 
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D n = nE n . (1.16) 



and 

Thus, (1.10) becomes 

V-D n = a n . (1.17) 
At oo, ^ gc ; = e = 0. Therefore k(oo) = 0, cr n (oo) = and 

D n (oo)=0. (1.18) 

Integrating (1.17) over all space we find, on account of (1.18), 

Ja n d N q = Je- 2S Xn d N q = 0. (1.19) 

Hence, the orthogonality condition (1.12) is now carried over to all n > 1. 

Expand e ~ s Xn in terms of the set of all eigenstates {ip a } of (1.1). We have, 
because of (1.19), 

e~ S Xn = E c a {n)ip a (1.20) 

with 

c «( n ) = 1? E^ e « ' c «( n ~ X )- ( L21 ) 

For problems like the one-dimensional quartic potential (1.5), when the cou- 
pling g 2 is large, only the first excited state a = 1 has an excitation energy 
e = Ei — E g d that is exponentially smaller than all other E a — E g ^. Thus, 
the corresponding c\(n) becomes exponentially large compared with all other 
c a ^i(n), in accordance with (1.21). Hence the iteration process (1.10) becomes 
rapidly convergent. 

For a large class of problems in which V(q) has high barriers and low valleys, 
there are often only a finite number of low-lying excited states. 

a = 1, 2, • - m (1.22) 

with excitation energies e a = E a — E g d comparable to each other, but expo- 
nentially smaller than those of a > m states. By maintaining the orthogonality 



relations between these m low-lying excited states, the iteration process (1.10) 
can be readily generalized to such problems. 

2. One Dimensional Problem 

In one dimension, we replace {qi\ by a single x. For simplicity, consider the 
special case 

V(x) = V(-x). (2.1) 

Thus, the groundstate wave function e~ s ^ is an even function of x and the 
first excited state e~ s ^xi x ) an °dd function. We assume that e~ s ^ is 
already known, e.g. by following the method discussed in Refs.[l-4]. We 
further assume that the first excited state has an excitation energy 

e«l, (2.2) 

as would be the case if V{x) is like the quartic potential (1.5) with large g. In 
one-dimension, (1.9) becomes 

The corresponding series of iterative equations (1.10) for n > 1 is 

1 ^ / -2S^Xn^ -2S /o A\ 

-2Tx {e ^ ) = £ » e Xn - 1 {2A) 

with xo a properly chosen trial function. As in (1.13), in order to separate Xn 
and e n from the ratio Xn/e n we choose a fixed point x° and set 

Xn(x°) = Xo(x°) (2.5) 

for all n. 

In terms of the electrostatic analog (1.14)-(1.16), write 

K = 6 , (T n = 6 n 6 Xn—1 

E n = -^x'n and D n = nE n . (2.6) 
with ' denoting d/dx. Thus, (1.17) and (1.18) become 

D' n (x) = a n (x) (2.7) 



£> n (±oo) = and a n (±oo) = 0. (2.8) 
From (2.7) and (2.8), we have 

oo 

D n (x) = -e„ / e- 2S W X n-i(z)dz, (2.9) 

x 

and therefore 

E n {x) = e 2S ^D n (x) (2.10) 

can also be readily expressed in terms of Xn-i- Since Xn{x) is odd in x, we 
need only to consider x > 0. From (2.6) and (2.10), we find 

X oo 

Xn(x) = 2e n j e 2S ^dy J e - 2S ^ X n-i(z)dz. (2.11) 
o y 

Since e~ s ^ is the groundstate, it has no zero at finite x. Thus, the factor 
e 2S(y) j n ^ e ^-integration of (2.11) is always finite. This is an important 
fact that enables us to extend the effectiveness of the iterative procedures of 
Refs.[l-5] for the groundstate to the low-lying excited state. 

3. An Analytically Soluble Example 

As an analytically soluble example, we consider the following simple one 
dimensional example 



V(x) = X5(x) + | 


o , 

oo , 


\x\ < 1 
\x\ > 1. 


(3.1) 


The unnormalized groundstate solution 


is 






e S(x) = < 


sinj9(l — x) 


sinj9(l + x) 


for 


< x < 1 
\x\ > 1 
- 1 < x < 


(3.2) 


with 


P = 7T 


-8 






and 








(3.3) 




A = —pcotp = 


(tt- 


5) cot 6. 





For very large A, p is very close to n, so 

6 < 1. (3.4) 
The exact lowest excited state is, for \x\ < 1 

e~ s x = sin7rx (3.5) 

with the excitation energy 

e = i(7r 2 -/) = ^-^ 2 «l, (3.6) 
and for x > 0, x ' s given by 

sinTrx 

= 7^ v 3.7) 

smp(l — x) 

In order to test the effectiveness of the iterative solution (2.11) and the sup- 
plementary condition (2.5), we shall start from the groundstate solution (3.2) 
and a trial function 

Xo = x. (3.8) 
Since (2.11) yields only the ratio Xn{x)/e n , we follow (2.5) by choosing 

x° = 1 (3.9) 

and therefore 

X„(l) = Xo(l) = 1. (3.10) 
It is clear that the trial function 

e~ s xo = x sinj9(l — x) (3.11) 

is not a very good guess of the first excited state (3.5), nor does (3.8) resembles 
(3.7), except that at x = 0, xo(0) = x(0) = 0. Nevertheless, we shall show 
that the iterative solution (2.11) with the supplementary condition (3.10) does 
lead to a rapidly convergent sequence. 

Substituting (3.2) and (3.8) into (2.11), we find, for x > 0, 

[- — sinp]xie~ S = sinpx — (smp)[— sinp(l — x) + x 2 cosj9(l — x)\. (3.12) 
-< l P 



The supplementary condition (3.10) leads to 



ei = — . 3.13) 

1 — J9COtj9 

For 5 = 7T — p small, (3.13) gives 

ei = 2tt(5 - 4£ 2 + 2(i + l)5 3 - 25 4 + 0(£ 5 ) (3.14) 

7T 3 

whereas the exact e is given by (3.6) with t\ = 2e. However, on second and 
third iterations, we find 

e 2 = 7T5 + (i - ^)<5 3 + (1 - 4,)<5 4 + 0(5 5 ) (3.15) 

Comparing to the exact value, the second order solution e 2 gives the correct 
tt5 and the third order 63 gives the correct first two terms: ttS — ^ 6 2 . 
For 6 = 0.1, (3.14)-(3.16) yield 

ei = 0.59086, e 2 = 0.31348, e 3 = 0.30924 (3.17) 

which can be compared with the exact solution 

e = 0.30916. 

The Xn(x) for n = 0, 1, 2, 3 are plotted in Fig.l, together with the exact 
solution x( x )- ' n f ac t< X2, X3 ar| d X are on the same curve. This shows the 
rapid convergence of the iterative process. 
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0.2 0.4 0.6 0.8 1 

Fig.l x ar| d Xn for n — 0, 1, 2, 3 for the soluble problem (Section 3). 
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4. One-dimensional Quartic Potential 



Consider the one-dimensional quartic potential 

V = 9 ^-lf (4.1) 

Assume that the groundstate wave function e has been obtained by using 
the method described in our recent paper[5]. The first excited state e~ s \ 
satisfies (2.3) and the corresponding n th order iterative solution is e ~ Xn, with 
Xn given by (2.11). When n = 0, we choose the trial function xo to be an 
odd function, with 

Xo(-x) = -Xo(x) 

and for x positive 

, x \ x(2 — x) for < x < 1 ,. n . 

»W = { 1 for x>l (42) 

As in (2.5) and (3.9)-(3.10), we choose x° = 1 and set 

X„(l) = Xo(l) = 1. (4.3) 

In the following, we also set 

g = 3. (4.4) 
Using the iteration method of Refs.[l,5], we find the groundstate energy to be 

E gd = 2.48291. (4.5) 

The lowest excitation energy e n for the first four iterations (n = 1, 2, 3, 4) 
based on (2.11) and (4.3) are 

ei = 0.41776, e 2 = 0.41367, 

e 3 = 0.413568, e 4 = 0.413568. (4.6) 

Thus, with three iterations, 63 is already accurate to seven significant figures. 
The corresponding wave functions Xn are plotted in Figure 2. The groundstate 
e~ s and the first excited state e ~ s x are given in Figure 3. To the same seven 
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significant figures the eigenvalue for the first excited state can be expressed as 

E dd = E gd + e 4 . 

To test the sensitivity to the choice x° we set, instead of (4.3), x° = 1/2 
and require 

Xn(i)=Xo(i). (4.7) 

The corresponding excitation energies e n of the lowest excited state for the 
first 4 iterations are, in place of (4.6), 

ei = 0.41363, e 2 = 0.41358 

^ = 0.413569. e.A = 0.413568. (4.8) 

Alt 
aft 
me 




0.5 1 1.5 2 2.5 3 

Fig. 2 Xn for = 0, 1, 2, 3, 4 for the quartic potential (4.1) with g — 3. 
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Appendix 



In our earlier work[6] an asymptotic expansion of the average E = ^{E g( i + 
E dd) and that of the difference A = \ e were obtained. To compare our 
present results with these asymptotic expansions, we take a larger coupling 
g = 8. This enables us to use the asymptotic expansion up to order 1/g 3 for 
the wave function (f>± of Ref.[6]. We find 



Using the expressions 



and 

of Ref.[6], we obtain 



E asy m P = 7.728854 . (Al) 



A = ? oo * - (A2) 
z / (f)\[x)dx 



A = - (A.3) 



tasymp — ^^asymp — 0.003027. (^-4) 

For g = 8, the iterative method of Ref.[5] gives 

E gd = 7.727340. (A5) 

The corresponding n th order iterative excitation energy e n of Sec. 4 for n = 
1 - 4 are 

ei = 0.00310125, e 2 = 0.00301796 

e 3 = 0.003017947, e 4 = 0.003017947. (A.6) 
Keeping the accuracy to seven significant figures, we find 

E = E gd + ^ e = 7.728849 (A7) 

and 

e = 0.003018. (A8) 

Thus, the asymptotic values E asymp and e asymp of (A.l) and (A. 4) compare fa- 
vorably well with E and e of (A. 7) and (A. 8). However, inclusions of still higher 
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X 



0.5 1 1.5 2 2.5 3 

Fig.4 Xn for n — 0, 1, 2, 3, 4 for the quartic potential (4.1) with g = 8. 

and higher order terms in the asymptotic expansion would lead to divergent 
results for E a symp and € aS y m p. 

The wave functions Xn{ x ) for g = 8 are plotted in Fig. 4. As we can see, 
the deviation of the Xn( x ) f rom 1 happens only at small x. 
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